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By using unitarized two-loop Chiral Perturbation Theory partial waves to describe pion-pion 
scattering we find that the dominant component of the lightest scalar meson does not follow the qq 
dependence on the number of colors that, in contrast, is obeyed by the lightest vectors. The method 
suggests that a subdominant qq component of the /o(600) possibly originates around 1 GeV. 

PACS numbers: 12.39.Fe,11.15.Pg,12.39.Mk,13.75.Lb,14.40.Cs 



so 

O 
O 

(N 
^— > 
O 

o 

o 



> 

OS 
CO 

o 

so 
O 

Oh! 
I 



X 



The lightest scalar mesons are a subject of a longstand- 
ing controversy that is recently receiving relevant contri- 
butions that could help settling the questions about their 
existence and nature. Experimentally, several analyses 
0, find poles for the Jo (600) and k, the lightest scalars 
with isospin and 1/2, respectively. The former is of in- 
terest for spectroscopy but also for understanding sponta- 
neous chiral symmetry breaking, since it has precisely the 
vacuum quantum numbers. On the theoretical side, the 
QCD chiral symmetry breaking pattern has been shown 
to lead to /n(600) and k poles in 7T7t and irK scattering 
[H IE LH IE Q • Concerning the spectroscopic classifica- 
tion, the caveat for most hadronic models is the difficulty 
to extract the quark and gluon composition without as- 
sumptions hard to justify within QCD. In contrast, when 
using fundamental degrees of freedom, i.e., in lattice or 
with QCD inspired potentials, other complications arise, 
related to chiral symmetry breaking, the use of actual 
quark masses or the physical pion or kaon masses and 
decay constants. All approaches are also complicated by 
the possible mixing of different states in the physical one. 
Most of these caveats are overcome in a recent approach 
7] (&] for a review) based on the pole dependence on the 
number of colors, N c , of meson-meson scattering within 
unitarized Chiral Perturbation Theory. 

The relevance of the large iV c expansion is that it 
provides an analytic approximation to QCD in the whole 
energy region and a clear identification of qq states, that 
become bound states as N c — > oo, and whose masses 
scale as O(l) and their widths as 0(1/N C ). Other kind 
of hadronic states may show different behaviors [10( . 

In order to avoid any spurious N c dependence in the 
hadronic description, we use Chiral Perturbation Theory 
(ChPT), which is the QCD low energy Effective Theory, 
and where the large iV c behavior is implemented system- 
atically. It is built as the most general derivative expan- 
sion of a Lagrangian in terms of tt, K and ry mesons 
compatible with the QCD symmetries, These particles 
are the Goldstone bosons associated to the spontaneous 
chiral symmetry breaking of massless QCD and are there- 
fore the lightest degrees of freedom. Actually, the u, d 
and s quark masses are non- vanishing but small enough 
to be treated as perturbations that give rise to 7r, K and 



r\ masses. Thus, ChPT is an expansion in powers of mo- 
menta and masses and, generically, its applicability is 
limited to a few hundred MeV above threshold. Each 
order is made of all possible terms multiplied by a "chi- 
ral" parameter. These Low Energy Constants (LECS) 
are renormalized to absorb loop divergences order by or- 
der, and once determined from experiment they can be 
used in any other pseudo-Goldstone boson amplitude. 

For our purposes, we are interested in meson- meson 
scattering amplitudes, since by unitarization they gener- 
ate dynamically resonances not initially present in ChPT 
U IE El Q El Indeed |E Ej, with the cou- 
pled channel Inverse Amplitude Method (I AM), the one 
loop ChPT meson-meson amplitudes describe data up 
to roughly yfs ~ 1.2 GeV and generate the p and K* 
vectors, as well as the /o(600), k, ao(980) and /o(980) 
scalars, and, most importantly, using LECS compatible 
with standard ChPT and therefore without any further 
assumption or source of spurious N c dependence. 

By scaling the one-loop ChPT parameters with their 
N c behavior, it was recently shown that the generated 
p and K* show the typical N c behavior of qq states, 
whereas the scalars are at odds with a dominant qq com- 
ponent. These results, confirmed by other methods [l5| . 
implied some cancellation between tree level diagrams 
proportional to LECS, and that loop diagrams with two 
intermediate mesons are very relevant in the generation 
of light scalars. But such loop diagrams are subdominant 
in the large N c counting and one could wonder about the 
stability under small changes in the LECS and about 
higher order ChPT corrections that could become larger 
than loop terms at sufficiently large N c , and reveal the 
existence of subdominant qq components. 

Here we present a method to quantify the above state- 
ments, and generalize the approach of to two-loop 
ChPT and in particular to 7T7t scattering EEl- Despite 
the many second order parameters and their large uncer- 
tainties, the data can be well described and we find once 
more that the p(770) behaves as qq with N c , whereas the 
Jo (600) main component does not behave as such. Fur- 
thermore, with the second order calculation a dominant 
qq behavior cannot be imposed on the /o(600) and the 
p simultaneously, but a subdominant qq Jo (600) compo- 
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TABLE I: Sample of LECS central values. The fifth and ninth 
columns, whose uncertainties roughly represent the previous 
sample, are used in our fits in the text to calculate xiscs- 



nent seems to arise at larger iV c around 1 GeV. 

Thus, at leading order, the only parameter is the pion 
decay constant in the chiral limit, /q = 0(\fW^) : fixed by 
the spontaneous symmetry breaking scale 47r/o — 1 GeV. 
Indeed, ChPT 7T7T scattering amplitudes are expanded as 
t ~ t 2 + t 4 + t 6 + ... with t k = O((p/ATTf ) k ) and are, 
generically, 0(1/N C ). In particular, the LECS appearing 
in 7T7T scattering at 0(p 4 ) [llj, all scale as 0(N C ). For 
simplicity we use the SU{2) notation, l?, since we are 
only dealing with 7T7r scattering (see the last reference in 
11] for a translation to SU(3)). In Table I we give a sam- 
ple of l\ sets given in the literature, whose differences we 
take as systematic uncertainties for the set we use in our 
fits below. In Table II we also list the six 0(p 6 ) constants 
that appear in 7T7t scattering, denoted rj. They all count 
as 0(N%). Those values are just estimates assuming they 
are saturated by the multiplets of the lightest (predomi- 
nantly vector) resonances. This hypothesis works well at 
0(p ) |2fJ, but for 0(p 6 ) is probably just correct within 
an order of magnitude |l6( and we conservatively assign 
a 100% uncertainty. 

The large N c counting does not specify at what renor- 
malization scale p, it applies, thus becoming an uncer- 
tainty studied in Q for the one-loop LECS. For the fj, 
the scale dependence is much more cumbersome and has 
not been written explicitly. Nevertheless, in both cases 
it is subleading in 1/N C1 and given the fact that we have 
a 100% error on the should be well within errors for 
our fits. Hence, we do not perform such analysis here, 
simply setting p — 770 MeV, as usual. 

Next, resonances can be found as poles in partial wave 
amplitudes tu of isospin I and angular momentum J 
that, in the elastic regime satisfy the unitarity condition: 

lmt = o-\t\ 2 ^lm- = -a^t = — — (1) 

t H,et~ L — ia 

where a is the known two-meson phase space and we have 
omitted the 13 indices for brevity. Note that ChPT ex- 
pansions violate exact unitarity, since in the first Eq.JU, 
the highest power of momenta on the right hand is twice 
that on the left. Unitarity is only satisfied perturbatively 

lmi 2 = 0, lmt 4 = atl, Imt 6 = a2t 2 Ret4, . . . (2) 



If we replace in Eq.QJ Ret -1 by its ChPT approxima- 
tion we get the Inverse Amplitude Method (IAM), that 
satisfies elastic unitarity exactly. At 0(p 4 ) it reads, 

t~t 2 2 /{t 2 -U), (3) 

and its fit to "data only" is listed in Table II, in the 
SU(2) notation. The fit quality is remarkable, given the 
huge systematic uncertainties, (conservatively ±5° and 
5% error for the Jo (600) channel) and we refer to H 
[hi ] for details and figures with a comparison with data. 
Using Eqs.JU and ©, the 0(p 6 ) IAM 0,113 reads 

t~tl/{t 2 -t A + tl/t 2 -t 6 ), (4) 

that recovers the 0(p 6 ) ChPT expansion at low ener- 
gies and describes well elastic 7T7t scattering data |2*lj . 
In addition, the IAM has a right cut that defines two 
Riemann sheets. In the second sheet we find poles asso- 
ciated to resonances; in particular, for the p(770) in the 
(I, J) = (1,1) channel and for the / (600) in the (0,0) 
one. For narrow resonances, T << M, the pole position 
is related to its mass and width as ^/s po i c ~ M — iT/2, 
and we keep this as a definition for the wide /o(600), 
whose M ~ 400 - 500 MeV and T - 400 - 600 MeV. 

By scaling the previous parameters with their dom- 
inant N c behavior, namely, /jy c - * fy/N c /3, l\ N — > 
l\N c /3 and r it N c — > ri(N c /3) 2 , we obtain the large N c 
dependence of M^ c and Tjv c of the p and /o(600) poles 
generated by the IAM. If a resonance is predominantly a 
cm, Mm ~ 0(1) and Tn ~ 0(1/N C ), and so it was shown 
|3 that the 0(p 4 ) IAM reproduced remarkably well that 
behavior for the p(770) and K * (892), two well established 
qq mesons. This is the expected behavior if in Eq.Q one 
neglects the two-meson loop terms, which are subleading 
at large N c with respect to 0(p 4 ) LECS contributions. 

In contrast, the lightest scalars follow a qualitatively 
different behavior. Loop diagrams, instead of the 0(p 4 ) 
LECS terms, play a relevant role in determining the 
scalar pole position. This is nothing but the well known 
fact that light scalars are dynamically generated by 
the resummation in Eq.@ of two-meson loop diagrams 
II II II 111 However, although relevant at N c = 3, 
loop diagrams are suppressed by l/N c compared to tree 
level terms with LECS, and the 0(p 6 ) terms could be- 
come bigger at some larger N c , where the 0(p 4 ) N c re- 
sults should no longer be trusted. For that reason it is 
important to check the 0(p 6 ) IAM: it should give small 
corrections to the 0(p 4 ) close to N c — 3, but it may de- 
viate at larger N c and even unveil some subdominant qq 
component. 

Still, before scaling the 0(p 6 ) IAM, let us first note 
that Mm c = O(l) and rjv e = 0(1/N C ) is only the leading 
qq scaling. Taking into account subleading uncertainties, 
to consider a resonance a qq state, it is enough that 

=m(i + —\ rf? = — ( 1 + — V (5) 

^ V AW Wc N c \ N C J' W 
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were M and T are unknown but N c independent and the 
subleading terms have been gathered in cm, er, which are 
0(1). Thus, for a qq state, the expected M Na and rjv e can 
be obtained from those at N c — 1 generated by the IAM, 
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Note the qq index for all quantities obtained assuming a 
qq behavior. We refer the values at N c to those at N c — 1 
because then we will be able to calculate from what N c 
value a resonance starts behaving as a qq, which is of 
interest in order to look for subdominant qq components. 
Thus, we can now define an averaged x| g to measure how 
close a resonance is to a qq behavior, using as uncertainty 
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Since at N c = 3 we expect generically 30% uncertain- 
ties we take cm = er = 1- Let us note that AM, and 
even faster AT, tend to zero for large N c and eventu- 
ally become smaller than our precision determining the 
pole position, 1 MeV, which we add as a systematic er- 
ror. When a state is predominantly qq, it should follow 
Eq.® and Xq q ^ 1- Otherwise Xq q ^ 1- Note that n 
should not be too far from 3, since we are looking for the 
N c behavior of the physical state. If we took n too large, 
we could be changing radically the original mixture of 
the observed state and for sufficiently large N c even the 
tiniest qq component could become dominant over the 
rest j2^. Therefore, our method first determines the be- 
havior of the resonance dominant component, but also, 
when qq is not dominant, the N c at which it becomes so. 

Furthermore, by minimizing its x? we can constrain 



a state to follow the qq behavior. 
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we will mini- 



mize Xdata + xlscs P 1us the Xqq ° f either of the p, or 
the /o(600), or both. The averaged Xlecs measures how 
far the fitted LECS are from their typical values in the 
tables and stabilizes them. Note that tttt scattering data 
are poor, with large systematic uncertainties and not very 
sensitive to some of the individual parameters but to their 
combinations, thus producing large correlations, driving 
the LECS away from the typical values for tiny improve- 
ments in the x 2 ; particularly at 0(p 6 ). We will provide 
the Xdatw xIecs> and Xq q , divided by the number of 
data points, the number of LECS and 2n, respectively. 

Thus, in Table II we show three 0(p 4 ) fits: to data 
only, constrained to a p qq hypothesis, or constraining 
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FABLE II: IAM fits to data or constrained to a qq N c behav- 
ior for the p and /o(600). In boldface the \ 2 minimized on 
each fit. For the 0(p 6 ) fits we also provide the ri estimates 

El 



the /o(600) to be a qq. We list for each fit the different 
X 2 described above and we see that our approach clearly 
identifies the p as a qq, since x 2 q q — 0.25. In contrast, 
Xq q > 125 for the /o(600), even if we constrain the fit to 
minimize also Xq q f° r the /o(600) (at the price of a higher 
x\ecs = 5.6). This is the quantitative statement of the 
0(p 4 ) results in 0, @ where it was concluded that the 
main component of the /o(600) was not qq. 

Unfortunately, the 0(p 6 ) analysis has a large freedom 
and thus Xlecs pl & y s a relevant role to stabilize the 
fit, but keeping in mind that the uncertainties were 
arbitrarily chosen to be 100%. In Table II and Fig. 1 we 
show three 0(p 6 ) fits: constraining the p as a qq (Fig. 1. 
Top), or the / (600) (Fig. 1. Center) or both (Fig. 1, 
Bottom). As expected, the 0(p 6 ) results are consistent 
with those at 0(p 4 ) not far from N c = 3 0, but for 
the scalar channel they deviate around N c ~ 8. 

In particular, in the "p as qq" fit a qq dominant nature 
comes out neatly for the p, whose x? ~ 0.9, but is dis- 
carded for the /o(600), since its Xq q — 15 and Fig.l shows 
that its mass and width both rise when N c increases not 
too far from real life, N c = 3. However, for N c > 8 the 
mass tends to a constant around 1 GeV and the width 
decreases, but not with a l/N c scaling. This suggests 
a mixing with a qq subdominant component, arising as 
loop-diagrams become more suppressed at large N c . 

One might wonder if the /o(600) could also be forced to 
behave predominantly as a qq. Thus we made a "/o(600) 
as a qq" constrained fit (Fig. 1, Center). The price to 
pay is a deterioration of the Xdata an d an unacceptable p 
qq behavior, since its x? ~ 2. Still, the /o(600) Xq q de- 
creases only to 3.5 (for 34 N c points). This extreme case 
allows us to conclude that the 0(p 6 ) calculation cannot 
accommodate a qq dominant component for the fo (600) . 
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FIG. 1: Mass and width N c behavior of the p and /o(600) 
from an 0(p 6 ) IAM data fit minimizing also the \qq- of the 
p (Top), of the /o(600) instead of the p (Center) of both the 
/o(600) and p (Bottom). 



Finally, we have studied how much of a subdominant qq 
behavior the /o(600) can accommodate without spoiling 
that of the p. Hence, we have also minimized in the 
fit the \\q both for the p and /o(600) (Fig. 1, Bottom). 
The /o(600) still does not behave predominantly as a qq, 
since its Xq q — 4- However, it starts behaving as a qq, 
i.e., Xq q — 1) for JV C > 6. The qq behavior of the p only 
deteriorates a little, Xqq — l-^> an d should not be pushed 
much further. This result suggests that the subdominant 
mixing with a qq state around 1 GeV seen in the first fit, 
would become dominant around N c > 6, at best . 

In summary, we have presented a method to deter- 
mine quantitatively how close the N c dependence of a 
resonance pole is to a qq behavior. We have applied this 
measure to the poles generated in rnr scattering by unita- 
rized Chiral Perturbation Theory, which is the effective 
low energy theory of QCD and reproduces systematically 
its large N c expansion. The method is able to confirm the 
0(p 4 ) qualitative results 0>H|] 5 identifying the p as a qq 
state and showing that the /o(600) is at odds with a dom- 
inant qq component. We have extended the method to 
0(p e ) confirming the stability of our 0(p 4 ) conclusions, 
but also showing that a possible subdominant qq may 



originate around 1 GeV. This provides further support, 
based on the QCD N c dependence, to some models that 
generate the /o(600) from final state meson interactions, 
and locate a "preexisting" qq scalar nonet 0, 0| around 1 
GeV. The methods presented here should be easily gen- 
eralized to investigate the nature of other dynamically 
generated mesons |23l | and baryons [24| . 
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